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Wprowadzenie transformacji Liego do optyki ma niedługą historię sięgającą 
początku lat osiemdziesiątych [1]. Transformacje Liego są pomostem łączącym 
optykę oraz mechanikę analityczną, pozwalającym wykorzystać wiele dobrze opra- 
cowanych technik rachunkowych formalizmu hamiltonowskiego oraz metod te- 
orii grup. W artykule niniejszym (oraz następnych) przedstawimy opis pro- 
pagacji światła w ośrodkach o zmiennym współczynniku załamania, na pozio- 
mie optyki geometrycznej z uwzględnieniem aberracji [2] oraz specyfiki układów 
światłowodowych. W dalszej kolejności przewidywane są prace nad „kwantyzacja” 


hamiltonianu optycznego pozwalającą uwzględnić falową naturę światła. Hamilto- 
nowskie sformułowanie optyki geometrycznej opisuje promienie światła jako linie w 
optycznej przestrzeni fazowej p = p(z)i q. = q (z), gdzie odpowiednikiem 
czasu z mechaniki klasycznej jest współrzędna | z” mierzona wzdłuż osi optycznej 


układu. Współrzędne uogólnione qı i q2 wyznaczają punkt przecięcia się promienia 
światła w „chwili z” z płaszczyzną z = const. Pędy uogólnione pı i p2 otrzymuje 


się z wariacyjnej zasady Fermata 
r 
35=8 j ds =0, (1) 


* Praca finansowana z funduszy Problemu RR - I — 02 
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gdzie ds = ndl = a dz oznacza droge optyczna, a n(g,z) wspölczynnik 
cos E 

załamania światła w punkcie q . Kąt 8 oznacza tu kąt pomiędzy normalną do 

płaszczyzny (q1,q2) a kierunkiem promienia przecinajacego tę płaszczyznę. Po- 


a s dq; 
nieważ tangens kąta 6 można wyrazić przez „prędkości uogólnione” q; = GT według 


tg 6 = Y (41)? + (42)? (2) 


to funkcjonał (1) można przepisać w formie 
r 
S= j La, 4 z)dz, (3) 


gdzie 
L(a, 4,z)=n(q,z)v1+ (41)? + (92)? - (4) 
Z funkcji (4) otrzymujemy pędy uogólnione odpowiadające współrzędnym uogól- 
nionym qı , 42 “= j 
n z + 


AK ie ES 


Wzory (4) i (5) prowadzą bezpośrednio do hamiltonianu optycznego 


H=),p4-L=-Vn2-p, (6) 


gdzie p? = (pi)? + (p2)?. 

Zgodnie z ogólnymi zasadami formalizmu kanonicznego z warunku wariacyj- 
nego (1) otrzymujemy równania ruchu dla dowolnej funkcji f( p, q ) określonej 
na przestrzeni fazowej 


T = _(H,f)=-As, (7) 


gdzie H = {H,-} oznacza (operator) nawias Poissona. W ogólności rozwiązanie 
równania (7) z hamiltonianem (6) dla dowolnie zmieniającego się współczynnika 
załamania n( q ,z) nie jest możliwe. Należy więc stosować rozwinięcia hamilto- 
nianu (6) w szereg do określonego rzędu. 

Doskonałym narzędziem pozwalającym otrzymywać przybliżone rozwiązania 
równania typu (7) dla złożonych układów optycznych jest analiza beotidgguno wa. 
Rozważmy nieskończoną algebre jednomianów dowolnego stopnia ZN = pę*pź pa De 
N = m + ną + m, + ma, z działaniem Liego w postaci nawiasu Poissona [2]. 
algebrą tą można związać algebrę operatorów działających w przestrzeni AS 
Operator odpowiadający funkcji f( p, q ) definiujemy wzorem 


Mo, ati ah 0 (8) 


Korzystając z tożsamości Jacobiego można pokazać, że zachodzi ważna równość 


[Aal = (4,9) 7 . (9) « 
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gdzie [7,5] = f3-— gf oznacza komutator, a { , }” operator ((f,9),'). 
Bezpośrednio z definicji nawiasów Poissona można łatwo sprawdzić, ze poniższy 
nawias Poissona 


{zYı, ZN3) = ZN + ZU + Z0 4 ZU; N=M+N-2, (10) 


jest wielomianem złożonym z jednomianów stopnia N = N + MNM -2 W 
szczególności dla N = N3 = 2 otrzymujemy dziesięciowymiarową algebrę 
generowaną przez jednomiany: 


PiPjy Wdy s Pili (i,j =1, 2) , (11) 


którą można zidentyfikować z algebrą Liego grupy symplektycznej Sp(4, R). Łatwo 
można sprawdzić, że jednomiany 


xi = Pp?=(p)’+(p)’, x605 P-a =maito, x= q*= KA te 

12 
także tworzą układ zamknięty względem nawiasu Poissona generując podalgebrę 
sp(2, R). Na dodatek, jednomian 


ma = q1P2 — 92P1 (13) 


generujący grupę SO(2) komutuje z jednomianami (12). Własność ta pozwala 
napisać następujący łańcuch grupowy 


Sp(4, R) > Sp(2, R) x SO(2) (14) 


odpowiedni do opisu układów optycznych o symetrii osiowej. Wynika to z faktu, 
że grupa SO(2) generowana przez operator ma = (ms, -) jest grupą obrotów prze- 
strzeni fazowej wokół osi optycznej „z” i dla hamiltonianu (6) niezmienniczego 
względem transformacji SO(2) układ optyczny musi posiadać oś symetrii. Jedno- 
cześnie należy zauważyć, że w tym przypadku hamiltonian optyczny (6) jest funkcją 
tylko elementów algebry sp(2, R). 

Uwzględnienie wyższych rzędów rozwinięcia hamiltonianu (6) niż drugi wy- 
maga rozszerzenia algebry sp(4, R) (dla układów osiowo symetrycznych sp(2, R)) 
o algebrę jednomianów wyższych stopni. Jak widać z zależności (10) jakiekolwiek 
wyjście poza algebrę sp(4, R) wymaga włączenia do rachunku jednomianów dowol- 
nego stopnia. Problem ten można rozwiązać tworząc algebrę ilorazową. Należy 
zauważyć, ze dla dowolnego ustalonego N, > 2 jednomiany stopnia wyższego niż 
Na + 1 tworzą nieskończoną algebrę, którą oznaczymy a>N<+!, Algebra ta stanowi 
ideał pełnej algebry jednomianów a”, a więc można utworzyć algebrę ilorazową 
opisującą nieliniowe efekty optyczne [2] 


aNati — a> JaN«t1 - a (15) 


Mówiąc inaczej, algebra aN« jest utworzona ze wszystkich jednomianów Z?, Z°,..., 


ZNe+ z iloczynem Liego przyjmowanym jako zero, gdy otrzymane jednomiany są 
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stopnia wyższego niż N, +1. Dla układów osiowo symetrycznych algebra (15) ulega 
zawężeniu tylko do jednomianów złożonych z elementów algebry sp(2, R). 

Zasadniczą ideą dalszego postępowania jest użycie przybliżenia hamiltonianu (6) 
wyznaczonego przez algebrę aN*. W tym przybliżeniu hamiltonian H staje sie ele- 
mentem algebry aNe, skąd wynika, że ewolucję promienia świetlnego można opisać 
przy pomocy transformacji należących do grupy Liego AY» odpowiadającej algebrze 
aN«, Niech f( p, q ;z) oznacza interesującą nas funkcję określoną na optycznej 
przestrzeni fazowej. Z teorii równań kanonicznych wiadomo, że dla równania (7) 
istnieje operator ewolucji Gy(z) zdefiniowany wzorem 


f(p,a;0)*/f(P, a;2)=Gul)f(p, 450) = f( Pp’, a';0) (16) 


Wstawiajac definicje (16) do röwnania (7) otrzymujemy równanie na operator ewo- 
lucji promienia świetlnego 


Cute) = —AGute). (17) 


Z drugiej strony, dla każdego 2, Gy(z) musi być elementem grupy transforma- 
cji optycznych AN*. Wymóg ten oznacza istnienie właściwej dla rozwiązywanego 
problemu trajektorii a(z) na rozmaitości grupowej AŃ* takiej, że 


AN* 3 G(a(z)) = Gp(z). (18) 


Wstawiając równość (18) do równania (17) otrzymujemy równanie określające 
trjektorię, wzdłuż której ewoluuje promień świetlny 


d > 
q Hala) = -Hgla(z)) , (19) 


tzn. operator A jest wektorem stycznym do linii a = a(z). W najprostszej a za- 
razem najważniejszej praktycznie sytuacji, gdy światłowód jest jednorodny wzdłuż 
osi optycznej „z”, hamiltonian optyczny jest niezależny od z i równanie (19) można 
bezpośrednio scałkować otrzymując eksponencjalny operator ewolucji promienia 
świetlnego 


Gy (z) = G(a(z)) = e77". (20) 
Z (20) wynika, iż w tym przypadku operator ewolucji tworzy jednoparametrową 
podgrupę transformacji optycznych, którą można skonstruować przy pomocy stan- 
dardowych metod teorii grup. W ogólnym przypadku hamiltonianu optycznego 
zależnego od z należy rozwiązać równanie (19) [2]. 

Dużą zaletą opisanej powyżej metody jest jej niezależność od warunków początko- 
wych, pozwalająca rozważać zarówno światłowody długie jak i krótkie oraz systemy 
świutłowodów połączonych w różny sposób. Wynika to z faktu, 12 można znaleźć 
operatory ewolucji promienia świetlnego (lub jakiejkolwiek funkcji opisującej wiązki 
promieni) dla każdego podukładu i w końcowym etapie rachunku je poskładać. 
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Jako przykład ilustrujący opisaną wyżej metodę postępowania, będący zara- 
zem podstawą dalszych rozważań, przebadajmy osiowo symetryczny światłowód 
gradientowy o współczynniku załamania określonym funkcją 


n(q) = no—vq? (= q7, no,v= const ). (21) 


W celu uproszczenia rachunków przyjmijmy przybliżenie optyki gaussowskiej, dla 
której grupą transformacji optycznych jest grupa symplektyczna Sp(2, R). W tym 
przybliżeniu hamiltonian optyczny zapisujemy w postaci 


p? 
H = — + vg? 3 (22) 
2no 
gdyż generatorami grupy Sp(2, R) sa trzy jednomiany (12), a hamiltonian optyczny 
musi być ich kombinacją liniową. Ogólną postać elementów grupy Sp(2, R), jako 
grupy transformacji optycznych w przestrzeni fazowej można zapisać w postaci 
wykładniczej [2] 


G(M( u )) = expl(u-1£1, +uoXlı + ugt), (23) 


gdzie M(u) oznacza podstawową, dwuwymiarową reprezentację macierzową 
grupy Sp(2, R) sparametryzowaną zmiennymi um; m = —1,0,1. Aby znaleźć po- 
stać M (u) wyznaczmy działanie generatorów Qi, na punkty przestrzeni fazowej 


ale port © = [3 | 
|; |-| 4 |-|o 4] q 
AAA O O 04 p 
E > EN a 
R" o TS 
sri blac W. a oromim i | 


Zgodnie z ogólną definicją działania operatorów grupy transformacji na przestrze- 
niach funkcji [3] macierzowymi odpowiednikami generatorów X}, są 


zj | 
al A 
ooo] (25) 
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Stąd wynika, że podstawowa, dwuwymiarowa reprezentacja grupy Sp(2, R). może 
być zapisana w formie 


LG 27 | il 03 
M(u) =exp (ua | | 0 yw 


ee Nuke 


2u1 uo 
Korzystając ze znanych wzorów [4] macierz (26) można przepisać w prostszej do 
zastosowań postaci 


(26) 


p sh w sh w 
ch w — ug zu 
M( u ) = sh ae sh w , (27) 
2u1 ch w + u —— 


gdzie w = +y uĝ — 4uju_1. 
Wynik ten jest identyczny z rezultatem podanym w pracy [2]. 

Po tym przygotowaniu możemy podać macierzową reprezentację hamiltonianu 
(22), który w postaci operatorowej zapisujemy jako kombinację liniową generatorów 
Im 

Ls 3 
A= Ino? + ux). (28) 
Korzystajac ze wzorów (25) otrzymujemy 


0 —2y 
A —H= RA 0 € sp(2, R) (29) 
no | 
Z równania (20) otrzymujemy, że operator ewolucji układu optycznego opisywanego 
równaniem (7) jest jednoparametrową podgrupą grupy transformacji optycznych 
generowaną przez hamiltonian (29) 


Gy(z) = Ë. (30) 


Analogicznie jak w przypadku macierzy (26) macierz (30) można zapisać w 
zwartej formie 


COS KZ —KNo sin KZ 
= R 
Gu(2) — sin Kz COS KZ d (31) 
KNO ‘ 


gdzie x = + = Stąd natychmiast otrzymujemy równanie opisujące dynamikę 
no . 


promienia świetlnego wchodzącego do światłowodu w punkcie q o pędzie (kie- 
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runku) p 
i 7 - " p cosK2— q Knosin xz j 
A E pee (82) 
q (a) | E q O 


Wzór (32) opisuje oscylacyjny ruch promienia świetlnego w światłowodzie o 
współczynniku załamania określonym przez równość (21). Częstość oscylacji jest 
pierwiastkiem z podwojonego stosunku gradientu zmiany współczynnika załamania 
do jego wartości na osi światłowodu. W użytym przybliżeniu wzory (32) nie 
uwzgiędniają żadnych aberracji optycznych. Problem ten jest tematem następnych 
artykułów. 
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SUMMARY 


In the paper the principles of group theory analysis for light propagation in the optical wave- 
guides with a graded refraction coefficient is presented. 
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